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ABSTRACT 

Photon  correlations  and  statistics  of  phase -conjugated  resonance 
fluorescence  of  a  two-state  atom  is  considered.  The  Q-factor,  as  a  function 
of  the  incident  laser  power,  the  deturning,  the  laser  linewidth  and  the  phase- 
conjugate  reflectivity,  has  been  calculated.  It  is  shown  that  for  small  and 
large  reflectivity  the  statistics  is  predominantly  sub-poissonian.  For  unit 
reflectivity  the  statistics  appears  to  be  exactly  poissom' an  (0  -  O'*  for  all 
values  of  the  optical  parameters. 


PACS:  42.65.H,  32.80,  52.50.D 
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1 .  Introduction 

Photon  counts,  as  measured  by  a  photomultiplier,  appear  as  random  events 
on  the  time  axis  FI, 2].  The  intensity  I(t)  is  defined  as  the  photon  counting 
rate  at  time  t,  and  therefore  I(t)dt  equals  the  average  number  of  detected 
photons  in  [  t ,  t  +  dt  ]  .  Then,  in  a  time  interval  [0,t]  the  average  number  of 
counts  is  given  by 

nit)  -  j  dt'  I(t')  .  (1) 

0 


For  stationary  radiation  the  intensity  is  time  independent,  and  we  have  /j(c)  = 

It.  The  statistical  fluctuations  in  the  number  of  counts  in  [O.t]  are 

represented  by  tlie  variance  a~  {t)  of  the  count  distribution.  In  the  case  of 

9 

pure  random  events  the  statistics  is  poissonian,  for  which  a  (c)  =  Ai(t),  and 
deviations  of  Poisson  statistics  are  most  conveniently  expressed  in  terms  of 
.Mandel's  Q-factor,  defined  as  [3] 


Q(t) 


(2) 


Obviously,  Q(t)  >  -  1,  and  for  Poisson  statistics  we  have  Q(t)  =■  0.  .4s 
pointed  out  by  Mandel  {3],  any  radiation  field  which  has  a  classical  analogue 
must  necessarily  have  a  positive  Q-factor  for  any  counting  interval  [0,t]. 
For  such  fields  >.hc  variance  is  larger  than  the  average,  and  the  statistics  is 
said  to  be  super  -  po  i  ss  on  ian  .  Conversely,  any  observati^..  of  a  negative  0- 
factor  would  indicate  the  essential  quantum  nature  of  the  detected  radiation. 


I 


L)  U 
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Such  sub - po i s s onian  statistics  was  predicted  [4]  and  demonstrated 
experimentally  [5]  in  single-atom  resonance  fluorescence. 

The  Q- factor  can  be  expressed  in  terms  of  the  two-photon  correlation 
function  I2(t^,t2).  By  definition,  I^Ct^, t2)dt^dt2  is  the  probability  for  the 
detection  of  a  photon  in  [t,  .t^+dt,  ]  and  a  photon  in  [  t- ,  t-+dt.,  ]  ,  irrespective 
of  photon  detections  at  other  times.  For  stationary  fields,  I2(t^,t2)  depends 
only  on  t^  and  t^  through  r  -  t^  -  c^,  and  the  Q-factor  assumes  the  form  [6] 

t 

Q(t)  -  ^  I  dr  ((t  -  r)l2(0.r)  -  rl^}  .  (3) 

0 

In  terms  of  the  normalized  correlation  function 

1.(0. r)  - 

h(T)  =  - 1 -  ,  (4) 

I 

the  Q-factor  can  be  written  as 

t 

Q(t)  =  ^  I  d^(t  -  r)h(r)  .  (5) 

0 

Of  particular  interest  is  the  limit  of  long  counting  times.  For  t  -►  «,  Eq. 
(5)  reduces  to 

CO 

Q(®)  -  21  I  dr  h(r)  ,  ("6) 

0 
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provided  that  the  limit  exists.  When  we  adopt  a  Laplace  transform 


03 

h(s)  -  I  dr  h(r)  ,  (7) 

0 


the  long-time  Q- factor  becomes 


Q(<x=)  -  2Ih(0)  .  (8) 


2 .  Phase-conjugated  fluorescence 

We  consider  a  two-state  atom,  with  excited  state  ]  e> ,  ground  state  |  g> 
and  level  separation  which  is  positioned  near  the  surface  of  a  four -wave 
mixing  phase  conjugator  (PC).  The  medium  is  pumped  by  two  strong 
counterpropagating  lasers  with  frequency  w.  A  laser  beam  with  frequency  is 
parallel  to  the  surface  and  irradiates  the  atom.  Resonance  fluorescence  is 
emitted  by  the  atom  and  detected  in  the  far  field  by  a  photomultiplier.  It 
can  be  shown  [7]  that  the  negative-frequency  part  of  the  detected  radiation  is 
proportional  to  the  Heisenberg  operator 

b(t)  -  d(t)  -  e^^‘^^Pd^(c)  ,  (9) 

where  d  =  |eXg|  is  the  atomic  raising  operator.  The  (complex)  number  P  is 
the  Fresnel  reflection  coefficient  for  a  plane  wave  with  frequency  ,  which 
is  incident  on  the  surface  of  the  PC  under  the  same  angle  as  the  location  of 
the  detector.  This  coefficient  P  depends  in  a  complicated  way  on  ,  the 
angle  cf  incidence  and  the  polarization  of  the  pump  beams  [8].  The  term  d(t) 
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in  Eq.  (9)  represents  ordinary  fluorescence,  and  the  term  proportional  to  P  is 
the  phase-conjugated  image  which  is  reflected  by  the  PC.  We  shall  assume  that 
the  medium  of  the  four-wave  mixer  is  transparent,  so  that  there  is  no  ordinary 
specular  reflection  of  the  incident  fluorescence  radiation.  In  terras  of  b(t), 
the  intensity  of  tht  detected  radiation  is 

I(t)  -  <|:<b(t)b(t)  *>:}.  ,  (10) 

and  the  two-photon  correlation  is  given  by 

12(^1. C2)  -  ?^<|:<b(t^)b(t2)b(t2)^b(tj^)^>:}-  .  (11) 

The  parameter  ^  contains  the  proportionality  factor  between  b(t)  and  the 
detected  field,  and  ^  is  proportional  to  the  detector  efficiency.  An  overall 
time  retardation  between  emission  and  detection  has  been  suppressed.  The 
finite  linewidth  of  the  driving  laser  will  be  considered  to  be  brought  about 
by  a  stochastically  fluctuating  phase.  The  notation  “j^.  . in  Eqs .  (10)  and 
(11)  indicates  an  average  over  the  random  laser  phase,  whereas  the  notation 
<.  .  .>  represents  a  quantum  average. 

Transforming  Eq.  (10)  to  the  Schrbdinger  picture  yields 

I(t)  -  ^<f:rrB(c)p(t):f  ,  (12) 

where  we  introduced  the  Liouville  operator  j^Ct),  which  is  defined  by  its 
action  on  an  arbitrary  Liouville  vector  II  according  to 
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]p(t)n  -  (d^  -  e’^^‘^'^P*d)n(d  -  (13) 

In  expression  (12),  p(t)  is  the  atomic  density  operator.  Its  time  evolution 
operator  will  be  indicated  by  the  Liouvillian  X(t,t'),  which  is  defined  as 

p(t)  -  X(t.c')p(c')  .  t  >  f  .  (14) 

With  this  notation,  the  two-photon  correlation  from  Eq .  (11)  becomes  in  the 

Schrodinger  picture 

=  ?V-"B(t2)X(t2,t^)p,(t^)p(t^):f,  ^^2  “  ^1  ' 

The  terms  inside  the  brackets  -f:.  .  in  both  Eq.  (12)  and  Eq .  (15)  are  randomly 
fluctuating  functions.  Only  their  average  over  the  stochastic  laser  phase 
will  reach  a  steady  state. 

3 .  Equation  of  motion 

The  Hamiltonian  of  the  two-state  atom  is  given  by 

H  =  Mcj  P  +  P  ,  (16) 

a  e  e  g  g 

in  terms  of  the  projectors  P  ■»  jeXej  and  P  =•  |gXg|  on  the  excited  state 

6  S 

and  ground  state,  respectively.  The  interaction  between  the  atom  and  the 
laser  field,  in  the  rotating- wave  approximation,  is  represented  by  the 


Hamiltonian 


7 


H  (c)  - 
ar 


-'2Wnc^' 


+  4>(t)) 


+  H.c. 


(17) 


where  H  is  the  (complex)  Rabi  frequency  of  the  dipole  coupling  and  4>(t)  is  the 
stochastic  laser  phase.  Then  the  equation  of  motion  for  the  atomic  density 
operator  p(t)  becomes 


"^dt 


[H. 


H  (c),p]  -  iy<rp 


(18) 


Spontaneous  decay  and  excitation  is  accounted  for  by  the  Liouvillian  f,  given 
by 


rn  =■  ~A  (P  n  +  np  -  2d^Id)  +  (p  n  +  np  -  2dJid^)  ,  (19) 

2eee  2ggg 

which  defines  its  action  on  an  arbitrary  Liouville  vector  H.  The  rate 

constants  A  and  .A  are  [9] 
e  g 

Ag  =■  A(1  +  i|p|-)  ,  Ag-JA|P|2  .  (20) 

with  A  the  Einstein  coefficient  for  spontaneous  decay  of  an  atom  in  empty 
space . 

Oscillations  with  the  laser  frequency  in  the  Hamiltonian  can  be 

eliminated  with  a  transformation.  With  the  Liouvillian  L  defined  as 

g 


L  n  -  [p  ,n] 

g  g 


(21) 


for  n  arbitrary,  the  transformed  density  operator  ait)  is  given  by  [10] 
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<7(t) 


- i(w  t  +  )L 
e  ^  ^  p(t) 


From  Eq .  (18)  we  then  obtain  the  equation  of  morion  for  a(t), 


(22) 


flfT  • 

i^  -  (L,  +  4>(t)L  -  ir)a  .  (23) 

dt  d  g 

Here,  the  dressed-atom  Liouvillian  equals 

L,n  =  AL  n  -  ^[nd  +  u*d^n]  ,  (24) 

d  g  2 

with  A  -  -  to  the  detuning  from  resonance. 

Equation  (23)  ic  a  stochastic  differential  equation  for  c7(t)  .  Ue  shall 
take  <^(t)  to  be  the  independent  -  increment  process  [llj,  which  has  the  phase- 
diffusion  model  as  its  gaussian  limit.  Then  Eq .  (23)  can  oe  solved  for  the 

average,  with  the  formal  result  [10] 

-i(L  -  iW  -  ir)(t  -  t  ) 

<j:a(t)|  =  e  °  -  (25) 

The  operator  W  accounts  for  the  phase  fluctuations  and  is  given  by 

W  -  ,  (26) 

g 

with  A  the  half-width  at  hal  f  -  maximum  of  the  Lorentzian  laser  profile.  The 
steady-state  value  •|:a(®):f  will  be  indicated  by  a,  and  is  the  solution  of 

(L^  -  iW  -  ir)o  -  0  .  (27) 

d 
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This  equation  is  easily  solved  for  the  matrix  elements  of  a.  For  the 
population  of  the  excited  state  we  obtain 

+  A|Pl^(A^+r7^) 

n  -  <e|CTle>  -  -  -  r  r  ,  (28) 

^  n  n  +  2A(A  +r;  ) 

o 

where  we  have  set  A  -  ^(A  +  A  )  ,  n  -  A  +  A  and  Q  -  IqI  .  The  population  of 

2  e  g  o'' 

the  ground  state  is  n  -  1  -  n  ,  arid  also  tne  coherence  <e  |  <7 1  g>  can  readily  be 

S  ® 

found  (but  is  not  needed  here) . 


d ,  Inters i tv 

With  Eq.  ''22),  expression  (12)  for  the  intensity  can  be  transformed  to 
the  c7  -  representation ,  and  with  Eq .  (13)  this  gives  four  separate  terms.  Then 
we  take  the  stochastic  average  and  the  limit  t  ®.  In  the  c- representation , 
the  terms  proportional  to  P  and  P*  acquire  phase  factors  of  the  form 
exp(±2i<p(t) )  .  With  the  identity  [12] 

lim  •|:e'’^'^‘''^^a(t) j  =  0  .  (29) 


these  cross  terms  vanish  identically  in  the  steady  state.  When  we  introduce 

the  Liouville  operators  and  ]R  as 

®  g 


d'nd  =  P  <a  n  e> 
g 


P  n  -  Ipi^dUd^ 

•  g  ' 


I^Pg<g|n|g> 


(30) 


p 


(31) 
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then  the  steady-state  intensity  can  be  written  as 


For  the  two  contributions  we  write 


(32) 


,  a  -  e,g  ,  (33) 

and  with  Eqs.  (30)  and  (31)  this  is 
The  total  intensity  is 


-  tn  and  I  -  F i P 
e  g 


respectively . 


The  significance  of  the  two  contributions  follows  from  Eqs.  (30)  and  (31). 
The  part  is  brought  about  by  the  action  of  on  the  density  operator  a, 

which  giv ‘t  "  '^glg^gl-  Therefore,  the  probability  for  the  emission  of 

an  "e-photon"  is  proportional  to  Che  population  of  the  excited  state,  and 

after  the  emission  the  atom  is  left  in  the  ground  si.'-te.  Similarly,  the 

action  of  IR  produces  a  "g-photon".  According  to  Eq .  (31),  the  probability 

g 

'^or  this  process  is  proportional  to  the  population  of  the  ground  state,  and 
after  Che  emission  the  atom  is  in  the  excited  state.  It  can  be  shown  [7]  that 
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this  stimulated  transition  actually  involves  a  three-photon  process;  twc 
photons  with  frequency  w  are  absorbed  and  a  photon  (the  g-photon)  with 
frequency  Iw  -  is  emitted  as  fluorescence.  An  atomic  transition  from  |  g> 
to  I  e>  then  guarantees  conservation  of  energy.  The  e-photons  are  ordinary 
fluorescence  in  an  j  e>  -*  |  g>  transition,  and  they  have  frequency  (in  the 
weak- field  limit). 

5 .  Two-photon  correlation 

The  two-photon  correlation  function  from  Eq .  (15)  can  be  worked  out  in 
the  same  way  as  the  intensity.  In  the  steady  state  (t^  -►  m,  r  =■  t^  -  t^ 
fixed)  the  cross  terms  vanish  identically,  and  we  obtain 

-i(L  -iW-iDr 

I^CO.r)  -  +  5^g)e  ^  •  (36) 

This  can  be  written  as 

a,/5 


in  terras  of  the  four  functions 


(r) 


^  -i(L  -iW-iDr 


a,P  -  e  or  g 


(38) 


From  the  interpretation  of  the 
the  intensity  of  /3-phorons  at  a 


R  operators  it  then  follows  that  equals 

time  r  after  the  detection  of  an  a-photon. 


From  Che  identity 
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-i(L^-iW-ir)r 

lim  e  n  -  aTrU 


for  arbitrary  n,  we  then  obtain 


lim  f^^(r)  -  ?rrJR  ^  -  I 

f -♦CO 


This  illustrates  that  for  a  long  delay  time  r  the  detection  rate  of  ;3-photons 

equals  the  uncorrelated  intensity  I..  With  Eq.  (37)  this  gives 

P 


lim  I  (O.r)  -  )  II  -  (I  +  I  -  I^ 

2  I  a  eg 


i.e.,  the  correlation  function  factorizes.  For  r  »  we  find 


f  ,(0)  -  f^^(O)  -  0  , 

QQ  gg 


f  (0)  =  ^  >  I 
eg  e 


f  (0)  -  ^|P|"  >  I 
ge  g 


Equation  (42)  expresses  that  the  probability  for  the  detection  of  an  e -photon 
immediately  after  the  detection  of  an  e-photon  is  zero  (antibunching).  This 
can  be  understood  from  the  fact  that  after  the  emission  of  an  e -photon  the 
atom  is  in  the  ground  state.  It  takes  a  finite  time  r  for  the  atom  to  make  a 
I g>  -»  j  e>  transition,  which  is  necessary  for  the  emission  of  a  subsequent  e- 


photon.  In  a  similar  way  it  follows  that  f  (0)  must  be  zero.  Equation  (43) 

gg 
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shows  chac  the  probability  for  the  detection  of  an  e-photon  immediately 

following  a  g-photon  is  larger  than  the  uncorrelated  probability  for  the 

detection  of  an  e-photon  (bunching).  This  follows  from  the  fact  that  after 

the  emission  of  an  e -photon  the  atom  is  in  the  ground  state,  rather  than  tr. 

This  enhances  the  probability  for  the  emission  of  a  g-photon,  which  is 

proportional  to  the  population  of  the  ground  state.  The  inequality  in  Eq. 

(44)  can  be  explained  in  a  similar  way.  From  Eqs .  (42) -(44)  and  with  n  +  n 

®  6 

+  1  we  find 


12(0, 0)  - 


(45) 


It  appears  that  l2(0,0)  does  not  depend  on  any  of  the  parameters  ,  A,  A  or 
i. 

The  Laplace  transform  of  Eq.  (38)  is 


”1  s  +  iL^  +  W  +  r  ’ 'a 

ct  a 


R  ;  , 


(46) 


in  terras  of  an  operator  inversion.  Working  out  this  expression  then  yields 
for  the  four  combinations 


(47) 


-  STiT  (4+A,)l<4  *  1)'  4^1 


(^8) 


^eg^"^  ’  lofry  ^  ^  (s+Ag)[(3  -H  7)^  +  A^j)  , 


(A9) 
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“  firir  +  n)  +  a^((s  +  +  a^d  , 

gg  sD(s)  2  o  e 


where 


D(s)  -  +  'l)  +  (s  +  2A)  [(s  +  fj)^  +  A^] 


With 


Eq.  (37)  we  can  then  construct  I^CO.s).  The  result,  however,  is  not  very 


transparent . 


6 .  Photon  statistics 


From  12(0,3)  and  with  Eq.  (4)  we  can  calculate  h(s)  .  Then  the  Q-factor 
follows  from  Eq .  (8),  with  result 


Q(«)  =  ^ 


n^(l+|P|^){(A-A)(rj^+A^)-4f)^A)+A|Pl^(3+lPl^)in^(r?^-A^)-(r,^+A^)^) 

{fi^/7+2A(f7^+A^))^ 

o 


4|p|^(r7W) 

'n^r,(l+lP|^)+AlP|^(r7^+A^)(3+iPl^-) 


For  |p1  -►0  this  reduces  to  the  Q-factor  of  a  free  atom  [6], 

Close  to  resonance  (A  -♦  0)  and  for  a  small  laser  linewidth  (X  -*  0)  the  Q- 
factor  can  be  written  as 


Q(«)  -  ^q(x.y)  , 


which  depends  basically  only  on  the  two  parameters 
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X  -  |P|^  .  y  -  n^/A^  .  (54) 

The  function  q(x,y)  is 

x(3  +  x){y  -  J(1  +  x)^}  -  |y(l  +  x)^  2x 

q(x,y)  -  - - - —  +  -  .  (55) 

{y  +  j  (1  +  x)^)^  y  +  ^  (3  +  X) 

The  sign  of  q(x,y)  then  determines  the  regions  of  sub-  and  super-poissonian 

statistics,  and  this  is  shown  in  Fig.  1.  It  can  be  checked  by  inspection  that 

q(x,y)  has  a  factor  1  -  x,  and  therefore  we  have  q(x,y)  -  0  for  x  -  1,  all  y. 

This  is  the  vertical  line  in  Fig.  1.  The  curve  in  Fig.  1  gives  the  second 

solution  of  q(x,y)  -  0.  On  the  y-axis  (PC  absent)  the  statistics  is  always 

3  1 

sub  -  po  is  sonian ,  and  q(0,y)  has  a  minimum  of  at  y  -  ~.  In  absence  of  the 

4 

laser  (x-axis) ,  q(x,0)  decreases  monotonically  from  “  at  x  -  0  to  -1  at  x  -♦  ®. 
The  function  q(x,y)  is  discontinuous  at  (x,y)  -  (0.0)  and  has  a  saddle  point 
at  (x,y)  =  (1,2). 

7 .  Unit  reflectivity 

2 

The  reflection  coefficient  |p|  is  proportional  to  the  square  of  the 

intensity  of  the  pump  lasers  of  the  four-wave  mixer,  and  can  therefore  have 

any  value.  In  particular,  reflectivities  larger  than  unity  have  been  obtained 

experimentally  [13-15].  An  interesting  special  case  is  unit  reflectivity,  for 
2 

which  jp]  -  1.  With  Eqs .  (47)  -  (51)  we  then  obtain 

ee  ge 


f 

eg 


(O  +  f 


gg 


(r) 


(57) 
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These  two  combinations  of  correlation  functions  turn  out  to  be  independent  of 
r.  The  intensity,  Eq.  (35),  reduces  to 

I  -  ^  ,  (58) 


and  the  two-photon  correlation  factors  as 


1^(0, r)  - 


for  all  r.  With  Eq.  (4)  this  yields  h(r)  -  0,  and  therefore 


(59) 


Q(t)  -  0 


(60) 


for  all  t  and  any  combination  of  parameters.  In  general,  it  is  not  necessary 
that  h(r)  is  identically  zero  in  order  for  the  statistics  to  be  poissonian; 


for  t  only  the  average 
Furthermore,  the  functions 
as  illustrated  in  Figs.  (2) 


of  h(r),  in  the  sense  of 

f  (r)  separately  have  a 

pa 

and  (3) . 


Eq.  (6),  has  to  vanish, 
non- trivial  t -dependence , 


8 .  Conclusions 

We  have  studied  the  photon  correlations  and  statistics  of  resonance 

fluorescence  radiation,  emitted  by  an  atom  near  the  surface  of  a  PC.  The  two- 

time  intensity  correlation  appeared  to  have  four  distinct  terms,  each  of  which 

is  proportional  to  a  function  f^  (t).  These  f  (r)'s  were  shown  to  have  the 

pet  pet 

significance  of  the  detection  rates  of  /3-photons  at  time  r.  after  the 
detection  of  an  a-photon  at  time  zero.  From  these  correlation  functions  we 
constructed  the  long-time  Q- factor,  as  given  by  Eq.  (52).  For  A  -  1  -  0,  Q(®) 


L7 


could  be  expressed  in  i  two-parameter  function  q(x,y) ,  and  Che  conditions  for 

sub-poissonian  statistics  (q  <  0)  were  represented  pictorially  in  Fig.  1,  For 

small  reflectivities  |P|  the  statistics  is  always  sub-poissonian;  for  |p|  — 

1  the  photons  have  Poisson  statistics  for  any  value  of  the  optical  parameters; 

2 

and  for  large  valuer  of  |p|  the  statistics  becomes  again  sub-poissonian.  It 

2 

follows  from  Eq.  (52)  that  for  |pj  -  «  the  Q- factor  reaches  its  ultimate 
lower  limit  of  Q('=°)  -  -f/A. 

Photon  correlations  and  statistics  of  phase -conjugated  resonance  fluorescence  should  be  amenable 
to  experimental  obsersation.  When  an  atomic  beam  and  a  laser  beam  are  projected  along  the  surface  of 
the  PC,  above  which  they  intersect  then  the  fluorescence  can  be  detected  in  a  direction  perpendicular  to 
the  surface.  A  complicadon,  however,  might  be  the  background  radiauon  which  is  emitted  spontaneously 
in  all  directions  by  the  nonlinear  medium. 


■Acknowledgmencs 

This  research  was  supported  by  che  Office  of  Naval  Research  and  the 
National  Science  Foundation  under  Grant  CHE-9016789. 


18 


References 

[1]  P.  L.  Kelley  and  U.  H.  Kleiner,  Phys .  Rev.  136  (1964)  A316 . 

[2]  R.  J.  Glauber,  Quantum  Optics  and  Electronics,  eds .  C.  DeWitt,  A. 
Blandin  and  C.  Cohen-Tannoudj i  (Gordon  and  Breach,  New  York,  1965). 


[3] 

L.  Mandel,  Opt.  Lett. 

4 

(1979)  205. 

(^] 

D.  Lenstra,  Phys.  Rev. 

,  A 

.  26  (1982)  3369. 

[5] 

R,  Short  and  L.  Mandel, 

Phys.  Rev.  Lett. 

51  (1983)  384. 

[6] 

H.  F. 

Arnoldus  and  G. 

Nienhuis,  Opt.  Acta  30  (1983)  1573. 

[7] 

H.  F. 

Arnoldus  and  T. 

F. 

George,  Phys.  Rev.  A  43  (1991) 

3675. 

[8] 

H.  F. 

Arnoldus  and  T. 

F. 

George,  J.  Mod. 

Opt.  36  (1989) 

31. 

[9] 

B.  H. 

W.  Hendriks  and 

G. 

Nienhuis,  Phys. 

Rev.  A  40  (1989)  1892. 

[10] 

H.  F. 

Arnoldus  and  G. 

Nienhuis,  J.  Phys 

.  B ;  At .  Mol . 

Phys .  16 

2325. 

[11]  N.  G.  van  Kampen,  Stochastic  Processes  in  Physics  and  Chemistry  (North- 
Holland,  Amsterdam,  1981),  p.  94. 

[12]  H.  F.  Arnoldus  and  T.  F.  George,  J.  Phys.  B;  At.  Mol.  Opt.  Phys.  (1991) 


in  press . 

[13]  D.  M.  Pepper,  D.  Fekete  and  A.  Yariv,  Appl.  Phys.  Lett.  33  (1978)  41. 

[14]  D.  M.  Bloom,  P.  F.  Liao  and  N.  P.  Economou,  Opt.  Lett.  2  (1978)  58. 


[15]  A.  Yariv,  IEEE  J.  Quant.  Electron.  QE  14  (1978)  650. 


19 


Figure  Captions 

Fig .  1  Regions  of  sub-poissonian  (-)  and  super-poisscnian  (+)  statistics  for 

A  -  A  -  0  as  a  function  of  the  phase -conjugate  reflectivity  x  and  the  relative 

laser  power  y.  The  indicated  signs  are  the  signs  of  the  function  q(x,y).  At 

(x,y)  -  (1,2),  q  has  a  saddle  point.  The  curved  line  approaches  the 

asymptotic  value  of  x  -  3  for  y  -►  ®,  and  q(x,y)  is  negative  for  all  x  >  3. 

Fig .  2  The  two  curves  represent  f  (t)/^  and  f  (t)/^  as  a  function  of  Ar  for 

ee  eg 

2  2  2 

A-A-O,  jp|  -1  and  -  2A  .  Hiese  parameters  correspond  to  the  saddle 
point  in  Fig.  1.  The  dashed  line  indicates  the  value  for  Ar  -♦  which  is 
for  both  curves.  The  overshoot  in  the  curve  of  f  (r)/^  below  the  asymptotic 
limit  is  a  remnant  of  Rabi  oscillations.  The  Rabi  frequency  for  these  values 
of  the  parameters  is  ^Jlk. 

Fig.  3  Functions  f^^(r)/^  and  £  (r)/|  for  the  same  parameters  as  in  Fig.  2, 

The  asymptotic  limit  here  is  as  shown  by  the  dashed  line. 
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